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An introduction to the CSL (Continuous Spontaneous Localization) theory of dynamical wave 
function collapse is provided, including a derivation of CSL from two postulates. There follows 
applications to a free particle, or to a 'small' rigid cluster of free particles, in a single wave-packet 
and in interfering packets. 
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I. INTRODUCTION 

Standard quantum theory is readily applied to mea- 
surement situations. This requires additional (ad hoc-for 
this case only) information be supplied for each situation. 
One may hope the theory could be extended to describe 
reality independent of experiment and without the need 
for ad hoc information. Given a state vector for any 
physical system, no matter how large or complex, one 
may hope for a theory that specifies the state vectors 
corresponding to the possible realizable states of nature, 
and their probabilities of realization. 

Why can't this be done with standard quantum the- 
ory? Given a state vector, 



I^M) = \ a n)(a n \ip,t), 



(1) 



one would like to say that the \a n ) correspond to the 
possible realizable states of nature, and that \(a n \tp,t)\ 2 
are their probabilities of realization. 

There are two problems with this. The first is the 
preferred basis problem. Why \a n ) and not another or- 
thonormal basis |6„)? No one has been able to specify 
the needed preferred basis. Of course, in experimental 
situations and special model situations this is possible, 
but it has not been possible in general. 

The second I like to call the hopping problem. If 
|(a ra |'!/',i)| 2 is the probability of \a n ) being realized in na- 
ture at time t, and \{a m \ip,t)\ 2 is the probability of \a m ) 
(m ^ n) being realized in nature at time t + dt, then 

\(a n \iP,t)\ 2 \(a m \iP,t + dt)\ 2 

is the probability of occurrence of both these events 
which, of course, is not seen. The difficulty is that quan- 
tum theory just gives the probabilities of these events: it 
does not give the probability of transitions between these 
events. 

The resolution chosen by the Founding Fathers was 
to restrict quantum theory to experimental situations, 
with empirically defined preferred bases, and to adopt 
the 'collapse postulate,' to ensure that a preferred basis 
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state, once chosen by nature, would remain chosen. They 
thereby relinquished the hope to describe reality. 

An alternative resolution, described here, does not give 
up on this hope. It alters quantum theory so that it 
may be realized. The Schrodinger equation is modified 
by adding a randomly fluctuating term to the Hamil- 
tonian, to account for the probabilistic behavior of na- 
ture. There is a preferred basis built into this term. The 
state vector dynamically collapses toward one of these 
basis states, very slowly for micro-objects, very rapidly 
for macro-objects. 

The preferred basis is essentially the mass density ba- 
sis. This choice could not be made for standard quan- 
tum theory, since collapse to position eigenstates gives 
the particles infinite energy. It works here because the 
collapse never goes all the way. First, it would take in- 
finite time for that to occur. Second, it doesn't happen 
because of the interaction of the usual Hamiltonian dy- 
namics with the collapse dynamics (e.g., the example in 
Section EES). 



II. DERIVING CSL 

For collapse dynamics, one wants the following behav- 
ior. 

Given the state vector Eq.([T]) at time t = 0, the squared 
amplitudes x n (t) = KanlV', i)| 2 should fluctuate until 
eventually one amplitude becomes equal to 1, and the 
rest equal to 0: that is the collapse. 

Moreover, with repeated evolutions, the nth amplitude 
should eventually reach 1 for a fraction x n (0) of the evo- 
lutions: that is the Born rule. 



A. Gambler's Ruin Game 

There is a rather precise analogy to this behavior, the 
'gambler's ruin game.'[l[ Consider two gamblers, with 
$100 between them. Gambler 1 starts with $Xi(0), gam- 
bler 2 with $A 2 (0) = $100- ATi(O). They toss a fair coin: 
'fair' is crucial, making the game what mathematicians 
call a 'Martingale.' Heads, gambler 1 gives a dollar to 
gambler 2, tails, the reverse. The amount each possesses, 
Xn(tk), fluctuates. Eventually, the game ends, with one 
gambler in possession of all the money. 
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Define P(Xi) as the conditional probability that gam- 
bler 1 eventually wins, given that he has $Xi. Then, 



P(Xi) = \P{X 1 - 1) 



1 



P(X X 



1) 



That is, when the coin is tossed, gambler 1 can either 
lose the toss but win from there, or win the toss and 
win from there. The solution of this difference equation 
is P(X{) = AXi + B (A and B are constants). With 
boundary conditions P(0) — and P(100) = 1, the solu- 
tion is P(Xi) = Xi/100. 

Define x n (tk) = X n (tk)/W0. Thus, we have seen that 
if gambler 1 starts out with a fraction xi(0) of the total 
amount of money, that is the fraction of repeated games 
he wins. 

This is just the behavior we want for a state vector 
which is the sum of two basis vectors \a±), {02) ■ Thus, 
the two basis vectors may be thought of as competing in 
a continuous version of the gambler's ruin game. 

For, the gambler's ruin game is a zero-sum, discrete, 
fair, random walk with absorbing barriers (at xi = 0, 1) 
in discrete time, for x n (tk) = X„(t/ C )/100. 

Collapse is a zero-sum, continuous, fair, random walk 
with absorbing barriers in continuous time, for x n (t) = 

IKhM)l 2 . 

The continuous limit of discrete random walk is Brow- 
nian motion. Thus, it is natural to regard the x n (t) as 
undergoing some kind of Brownian motion. We shall re- 
place the probabilities associated with a single coin toss, 
and with a sequence of coin tosses respectively by 



P(dB) = 
P(B(t)) = 



1 



yj2TT\dt 



1 



y/2TT\t 



-dB 2 /2\dt 



-B{t) 2 /2\t 



(2a) 
(2b) 



where A is a constant diffusion rate. Just as the gambler's 
ruin dollar count X n (tk) depends upon a sequence of coin 
tosses, so we shall take the collapse x n (t) to depend upon 
B(t). 



B. Postulates 

We shall derive the CSL collapse dynamics from two 
postulates: 

1) gambler's ruin behavior. 

2) a linear, real, Schrodinger equation. 

2) could use some explanatory remarks. 

In order to achieve the linearity part of 2), it shall be 
necessary to relinquish the condition that the state vector 
norm is 1. In standard quantum theory, the unit norm 
condition is mandatory, since it ensures that probabili- 
ties add up to 1. Here there is no need for the unit norm 
condition: probabilities are provided by P(t)'s probabil- 
ities. Physical information is carried by a state vector's 



direction in Hilbert space, not its norm. To calculate ex- 
pectation values, one can always normalize a state vector. 
The un- normalized state vector shall be denoted \<f>,t). 

Regarding the real part of 2), we wish to consider a dif- 
ferent dynamics from the usual Hamiiltonian dynamics 
d\ij},t) = —iHdt\ij),t). That utilizes a hermitian hamil- 
tonian H, so we shall consider dynamics of the form 
d\ijj,t) — H c dt\ijj,t) where the 'collapse hamiltonian' H c 
is hermitian (we note that the most general linear equa- 
tion is of the form d\ip,t) — —i(H + iH c )dt\ip,i)). 

Now, for the purposes of this derivation, for simplicity 
and appropriate to the generalization of a classical game, 
we shall further restrict the state vector components to 
be real numbers. Therefore, we must restrict H c to be a 
real symmetric operator in the chosen basis \a n ). Once 
we find the unique form of H c that allows postulate 1), 
it may readily be seen that that complex state vector 
components and the most general hermitian H c also allow 
postulate 1). Because the state vector components are 
real, we may write 



x„(t) 



(a n \4>,t) 2 



(a n \(j),t) i 



<<M<M> E, 



(3) 



Also in order to achieve the linearity part of 2), the 
Schrodinger equation shall be linear in a Brownian mo- 
tion B'(t). However, the relation assumed between B'{t) 
and B{t) is allowed to be non-linear in the state vector: 



dB'(t) = dB(t) + f{x)dt. 



(4) 



where / is an arbitrary real function of the x n (t). For, 
collapse violates the superposition principle, and this 
somehow requires a non-linearity. What 2) really means, 
then, is that we look to have a linear Schrodinger equa- 
tion and isolate all the non-linearity in the probability. 

The derivation, which is presented in the next few sec- 
tions, is rather lengthy, and it involves stochastic differ- 
ential equations-but the result does not. Some readers 
may wish to move immediately to section IIIII where the 
results obtained here are utilized. 

Postulate 1) is implemented by an ltd equation for 
x n (t): 



dx n (t) = b n (x)dB(t), 



(5) 



where b n is an arbitrary real function of the x m (ty&. 

We shall denote by an overline the ensemble average of 
a quantity, e.g., from Eq.([2a|l. dB = / d(dB)P(dB) = 0. 

It immediately follows from Eq.© that dx n (t) = 0. 
This says that it is a 'fair game' for each 'player' which, 
as we have said, is crucial for gambler's ruin behavior. [3j 
This is why Eq.© is chosen to be an Ito equation. 

The other necessary condition for gambler's ruin be- 
havior is the end-game condition, that 6„(x) = when 
one x n {t) is equal to 1 and the rest vanish. However, 
that does not need to be separately imposed since it au- 
tomatically occurs, as shall be seen. 
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Postulate 2) is implemented by the Stratonovich 
Schrodinger equation for the un-normalized state vector 
amplitudes: 

d\4>,t) = [RdB' + Sdt]\4>, t) 

where R and S are arbitrary symmetric real oper- 
ators, as discussed above. This is chosen to be a 
Stratonovich equation because the calculus manipula- 
tions (e.g., derivative of the product of functions) are 
the usual ones, which would not be the case for an Ito 
equation. 

However, a rather tedious calculation (Appendix A) 
shows that R and S have to be diagonal in the \a n ) basis 
if postulate 2) is to imply postulate 1). Therefore, we 
shall write this equation as 

d(a n \4>,t) = [a n dB' + f3 n dt](a n \(f>,t) (6) 
where a ni fi n are real constants. 



where c is an arbitrary constant. Putting Eq. (|10l) back 
into Eq.([9bf gives 

(fi n —fi-x) = ~c(a n -a-x)-\(a^-a 2 -x), or fi n = -ca n -Xa 2 L 

(11) 

the last step following from applying d/dxi to the first. 
Putting / from Eq.([T0|) and fi n from Eq.([ll]) into Eq.©, 
we find that the latter is independent of c. Therefore, c 
has no physical effect. Its only effect is to add a constant 
drift to B'{t) (see Eq.©), and since we are free to choose 
B'(t) to be as simple as possible, we may take c = 0. 

This concludes our derivation. We have found that 
postulates 1) and 2) lead to the Schrodinger equation © 
and the Probability Rule in the forms 

d{a n \(f>,t) = [a n dB' - Xa 2 n dt](a n \(t),t}, (12a) 
dB'(t) = dB(t)+2Xa-xdt, (12b) 

where the a n are completely arbitrary constants and x 
is given by Eq.([3]). 



C. Derivation 



D. Schrodinger Equation 



We proceed to find dx n (t) from Eq.Q. With use of 
Eqs.(j4]),([6]), we obtain the Stratonovich equation 

dx n (t) = 2{(a n — a ■ x.)dB 

+ [(a„ - a ■ x)/ + (fin - fi ■ x)]dt}x n (t). (7) 

where a ■ x = ct m x m (t) 

Now we may use the rule for converting a Stratonovich 
equation to an Ito equation, which in this case means 
adding 



Xdt 



E 2 « 



<\ m - a ■ x)- — 2{(a n - a ■ x) 

OXrn 



to the right side of Eq. (JTJ) . The result is 

dx n (t) — 2{(a n — a ■ x)dB 
+ [(a„ - a -x)f + (fi n -/3- x)]<ft 
+ [-a 2 • x + (a • x) 2 + (a„ - a • x) 2 ] Xdt}x n (t). (8) 

where a 2 -x = J2 m a m x m{t)- In order that the ltd Eq.® 
(consequence of postulate 2) agree with the Ito Eq.© 
(consequence of postulate 1), the coefficient of dB and 
the coefficient of dt in both equations must be equal, so 

6 n (x) = 2(a n - a ■ x)x„(t), (9a) 
(a n - a ■ x)/ + (fi n - fi ■ x) = 

X[a 2 ■ x - (a ■ x) 2 - (a n - a ■ x) 2 ] . (9b) 



To find / and fi n , operate on Eq. (19bl) with 
J2 n x n d 2 /dx 2 n . Remembering that ^2 n x n = 1, we ob- 
tain 



2a, 



d 



dx r 



■f = 4Aa 2 n , or / = 2Xa ■ x + c. (10) 



We may introduce an operator A defined by -A|a n ) = 
a n \a n ), in terms of which the Schrodinger equation (|12a[) 
may be written in basis-independent form, 

d\(f>,t) = [AdB'-XA 2 dt]\4>,t) or \4>,t) = e AB ' {t) - A2xt \(j>,0). 

(13) 

Since the \a n ) are eigenstates of A with eigenvalues a n , 
we could instead label them |a„). However, we shall 
rather replace a n by a n . 



E. Probability Rule 



Eq. (|12b[) is the same form whether an Ito or 
Stratonovich equation. It may be written as 



dB'(t) = dB(t)+2Xdt(A)(t) where (A)(t) 



(^t\A\^t) 



(4>,t\cf>,t) 

(14) 

We shall use it to find the expression for the probability 

ofB'(i). 

It follows from Eq.CEJ) that dB' {t) = 2Xdt(A)(t), 



(dB'(t)) 2 = Xdt, and all higher moments of dB' (t) can 
be neglected as they are of higher order than dt. There- 
fore, dB'(t) is a gaussian process, characterized by just 
these two moments. One may immediately check that 
the probability density, variously written as 



P\dB'(t)] 



1 jjgW U t \e 2AdB (t)-^xdtu A 

- g 2Xdt 



\/2irXdt 



1 IX *l (dB'(t)-2>.dtA) 2 

1 (4>,t\e 2A rft \fat) 



v / 2irXdt 



(15a) 



1 (dB'wf (4>,t+ dt\d)A + dt) , rlN 
e 2Xdt — — - , (15b) 



V2nXdt 



(ht\<i>,t) 
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has these two moments. 

Define N = t/dt. We multiply Ea.([i~5bl by itself N+ 1 
times with successively smaller values of t, obtaining the 
joint probability of the independent increments dB at 
successive values of t: 

N N 1 , 2 

P[dB'(t-ndt)] = Yl ■ - p — 2A * 



n=0 



n=0 



\Z2n\dt 



(0,O|0,O) 

Since (0, O|0, 0) = 1, and using Eq. (IT^|) to write 



.(16) 



M + di|0,t + dt) = (0,O|e 2AB '(*+ dt )- 2A2A (* +dt )|0,O) 
0\°l H e 2A ( s '(*+ dt -™ d *)- s '(*-" dt )- 2A2Adt )| (? !) ; 0) (17) 



(taking 5(0) = 0), it follows from Eas.([l6 l) .([iT | that the 
joint probability of the values of B at successive values 
of t is 

N 1 

P[B'{t + dt),B'{t),..,B'(dt)} = Yl -= 



„=o V^nXdl 

■<0, O|e- (B ' (t+t ""- dt) ^r" dt) ' 2Ad " )2 |0, 0).(18) 

We have written eLB'(i - ndt) = B'(t + dt- ndt) -B'(t- 
ndt). We have also written P as the joint probability of 
{B'{t + dt), B'{t), .., B'(dt)}, instead of the joint proba- 
bility of {dB'(t),dB'(t - dt), ...dB'(0)}, which we can do 
since the Jacobian determinant for the change of vari- 
ables is 1. 

Now, we arc interested in finding the probability 
P[B'(t + dt)}, regardless of what Brownian path leads 
to the value of B'(t + dt). To obtain this, we integrate 
Eq.(HHl) over all B n s except B'(t + dt). These integrals 
are easily done, since each B'(mdt) appears in just two 
(adjacent) gaussians in the product, and 



1 (B-l) 2 1 (B-a 2 f 

dB P 3oi p 2c 2 



\J2-kc\ 



V2TTC2 

1 



y/2ir( Cl + c 2 ) 
The result for P[B'(t + dt)} is 



(°i-°2r 



P 



v /2?rA(t + dt) 
1 

v /2?rA(t + dt) 



(B'(t + dt)-2A(t + rft)A)^ 

i, 0|e 2Ad* |0, 0) 



(B'(t + dt)) 2 , 

e jxd^- i + di|0 7 i + dtji9) 



(1191) . It is useful to incorporate the exponential factor in 
(|19p in the Schrodinger equation (and un- prime B, and 
replace t + dt by t) , so that its solution then becomes 

\^,t)=e-^ B ^ 2XtA ^\4>,Q), (20) 
leaving the Probability Rule in the simple form 



P[B{t)}dB(t) = 



dB(t) , 
\/2n\t 



,t\(j>,t). (21) 



Eas.(|20 l) ,(|2T |l define CSL. Everything that follows is 
based upon these two equations. 



A. Collapse mechanism 

Let's see how they contrive to collapse the state vector. 
Suppose the initial state vector is 



10,0) = ^c„|a„). 



(22) 



|0,O) is assumed normalized to 1, so X) n l c n| 2 = 1- 
Eas.pD]). (HU) then become 



|0, t) =5>„K)e-ro[ s ( 



t)-2\ta n ] 2 



(23) 



and 



P[B(t)]dB(t) = ^B=y2\c n \ 2 e-^ B ^- 2Xta ^ 2 . (24) 
V27rAt n 

Now, assume all the a n are unequal. Then, Eq. (|2~4")) 
describes a bunch of gaussians whose centers at 2Xta n 
drift further and further apart, while their widths yXt 
spread much more slowly. As t increases, the gaussians 
have less and less overlap. 

Then, for the set of £?(f)'s which lie within the mth 
gaussian, the state vector and probability arc, to an ex- 
cellent approximation (which becomes exact for t — > oo), 



\^,t)*c m \a m )e-^ B ^ 2Xt ^ 2 



(25) 



and 



P[B(t)]dB(t) 



dB(t) 
V2n\t 



Ic.fe-^W- 2 *' ™]". (26) 



For this set of i?(i)'s, the integrated probability in (|2"6")l 
is w |c m | 2 , giving the Born Rule. 

For any B(t) in this set, the normalized state vector is 



III. CSL 

We have arrived at CSL's two equations in the form of 
the Schrodinger equation (fT5)l and the Probability Rule 



U t ) = , fe£ « \a m ) 
giving the collapsed state vector . 



(27) 
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B. Refinements: Density Matrix 

The density matrix, as usual, is constructed from 
the state vectors and their associated probabilities, here 
given by Eqs.®,®: 



Pit) 



P[B(t)]dB(t) 



l<M(<M 
(4>,t\4>,t) 



dim u , u j<t>,t)(<t>,t 



— oo 



(<t>,t\cl>,t) 
|0,O)(0,O|e- 



dB ^ P -^iBM-2\tA] 2 U ^ U n | p - 3 l T [S(t)-2AtA] 2 



i[A L -A R y 



|^,0)(</»,0|. 



(28) 



In Eq. (|28| , Al,Ar mean that these operators act on the 
left or the right respectively of the initial density matrix 
|0,O)(0,O|. 

With the initial density matrix (|2"2")l , the matrix ele- 
ments of the density matrix at time t is found from 
to be 



(a n \p(t)\a m ) = c n c* m e " [o ™ am]2 



(29) 



showing how the off-diagonal elements decay while the 
diagonal elements remain constant, the collapse rate in- 
creasing as the eigenvalue differences increase. 

According to Eq. (|2"51) . the differential equation satisfied 
by the density matrix is 

^- = ~[A,[A,p(t)}}, (30) 
which is the simplest possible Lindblad equation. Q 



C. Refinements: Hamiltonian 

To add the Hamiltonian to the state vector dynamics, 
consider the evolution over an infinitesimal time interval: 
it and the probability rule become 

\^t) = e -iHdt-^ rt [dB{t)-2\dtA]*^ t _ dt) 

= e - d ^ H +^ w ^- 2XA ^\(t>,t-dt) (31a) 

(<M<M) dw(t) 



P(w)dw 



t>,t-dt\<j>,t-dt) y /2nX/dt 



(31b) 



where w(t) = dB(t)/dt is called white noise. 

Over a finite time interval, Eqs. (|31a]) . ()31b| imply 

\<f>,t) = re -tidt' {i H(t') + ^Mt')~2XAf} |^ 0) (32a) 



t-dt 



p{w)dw = (<m<m> n 



dw(t') 
•Jj, y/2-irX/dt 



(32b) 



where T is the time-ordering operator. 

To summarize, for each white noise function w(t) there 
is a corresponding state vector \4>,t) given by Eq. (|32al) . 



one of which is supposed to be realized in nature with 
probability (|32b|) . Generally, the hamiltonian evolution 
and the collapse-hamiltonian evolution compete against 
each other. This can give rise to effects which suggest ex- 
perimental tests of the collapse theory vis-a-vis standard 
quantum theory /collapse postulate. 

'White noise' was named after the sound which has all 
frequencies in equal amounts, in analogy to white light. 
It was named in a paper on the acoustics in airplanes, [f| 
where the authors wrote: 

That white noise is annoying needs little argument. No 
one has been found who really enjoys it. 

However, here it is enjoyed, in its role as the "chooser" 
of the collapsed state. 



D. Refinements: More Collapse-Generating 
Operators 

It is a straightforward generalization to describe col- 
lapse to a joint basis of operators A a which commute, 
[A a , A&] = 0. This requires one white noise function w a 
for each A a . The state vector evolution is 

|0, *> =Te~ Xo*'0*(0+AE> a (0-2AA«] 2 }|^ 0) (33) 
and the corresponding density matrix evolution is 

^ = -i[H,p(t)] [A a ,p(t)}}. (34) 

a 

The ensemble average of an operator O shall be denoted 
0(t) = TrOp(t), where Tr is the trace operation. Then, 
Eq-® gives 



dO{t) 
dt 



A 



-i[0,n](t)--J2[A a ,[A a ,0}}(t). (35) 



IV. NON-RELATIVISTIC CSL 

Finally, here is the CSL proposal to describe the 
non-relativistic world. @, The index a in Eq. ([33|) is 
changed to a continuum index x, so the 'chooser' u>(x, t), 
rather than being a set of random functions, is a random 
field: 

\<f>, t) = Te~ l So dt'H{t')--L. p o dt' J dx>(x',t')-2AA(x')] 2 |^ Q y 

(36)' 

The set of collapse-generating operators are mass-density 
operators, 'smeared' over a sphere of radius a: 

(37) 

Here, ££(z) is the creation operator for a particle of type 
n at z. m n is the mass of this particle and M is the 
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mass of a nucleon (say, the neutron). Thus, in ordinary 
matter, it is the nucleons which are mostly responsible for 
collapse. Experimental resultsQ have dictated that the 
effective collapse rate in Eg. ([5r?|) be mass-proportional, 
~ Xm n . 

Assuming the theory is correct, the parameter values of 
A, a should be determined by experiment. 8] Until then, 
we shall provisionally adopt the parameter values chosen 
by Ghirardi, Rimini and Weber|9J in their instantaneous 



collapse theory, A ps 10 _16 sec _1 ,a ~ 10 _5 cm. However, 
it should be mentioned that Adler[Io| has given an argu- 
ment for A to be as large as ps 10 _11 sec _1 . 

The density matrix evolution equation Q34p becomes, 
using (|571) . 

dp{t) , r M] A^7n fc m„ 

= -i[H,p{t)\ - - 



dt 



2 ^ M 2 

k.n 



(tto 2 ) 3 / 2 



A 
2 



Af 2 



dz / dz'e-^l z - z '] 2 [ei(z)a(z)[^(z')^(z'),PW]]-(38) 



V. FREE SMALL CLUMP 

In the rest of this paper, we shall illustrate CSL by dis- 
cussing the force-free behavior of the center of mass (cm) 
of a small (dimensions< a) clump of ordinary matter. 1 

For simplicity we shall neglect the electrons, and take 
there to be N nucleons, regarded as a single type of par- 
ticle, of mass m = M , which are very good approxima- 
tions for our calculations. Then, in the particle posi- 
tion basis, |x) = |xi, ....x,, ...x N ), using ££(z)£ fe (z)|x) = 
J2iLi — x 0l x )i Eq. (|38| becomes 



d(x|p(t)|x') 



dt 



= -i(x\[H,p(t)]W)-(x\P®W) 



N 



(39) 



Define the center of mass coordinate X = N 1 Yli=i x * 
and the relative coordinates y; = x, — X. Because it is a 



'small' clump, exp — [y— yj] 2 /4a 2 w 1 and exp — [X+y, — 
X' - yj] 2 /4a 2 re exp-[X - X'] 2 /4a 2 . With the density 
matrix assumed to have the form of the direct product 
of cm and internal coordinate density matrices, we can 
take the trace over the internal coordinates in Eq. ([39| to 
obtain the equation for the evolution of (X|p(t)|X'). 

It is useful to express this equation in operator form, 
writing the cm operator as X, and its conjugate cm mo- 
mentum as P: 



dp(t) 
dt 



P 2 








-XN 2 


l-e 



rV[X L -X H ] 



pit). 

(40) 

The associated state vector evolution equation (|36j) is 

|</> t) =re- i ^o d *'2OTP 3 -TC/o d *'/ dx '^( x '>*')-2Aiv^(x')] 3 | 9i 

(41) ' 

where 



A(x') ee 



1 



<'-Xl ; 



>a 2 ) 3 / 4 



(42) 



To illustrate the use of the collapse part of Eq. (|40|) . 
consider the initial wave function 



1 



|0,O) =^=[\L) 



\R) 



(43) 



where the states describe the clump to the left or right, 
with the two wave packet em's separated by the distance 
D. Then (ignoring the kinetic energy), 



d(L\p(t)\R) 
dt 



1 - e 



(L\p(t)\R). (44) 



Thus, for D >> a, the collapse is described by expo- 
nential decay of the off-diagonal density matrix element 
with characteristic time A -1 = 10 16 sec for a single nu- 
cleon. For a 10~ 5 cm cube of gold, where N ~ 10 8 , the 
characteristic collapse time is 1/XN 2 = lsec. 

One might very well extend this theory to include 
massless particles by replacing mass-density of A(x) in 
Eq. ([57j) by energy-density /c 2 . One might then regard it 
as holding in the co- moving frame of the universe [12|, or 
as the limit of a relativistic CSLILJ. 



VI. COLLAPSE OF A PACKET 

We shall consider how a single wave packet undergoes 
collapse. 



1 The internal excitation of the matter is not discussed here. Col- 
lapse narrows wave packets, resulting in atomic and nuclear 
'anomalous' excitation (i.e., collapse- generated, not predicted by 
standard quantum theory). 7]. Experimental limits on such ex- 
citation strongly suggests the effective mass-proportionality of 
the collapse rate, as we have mentioned. Incidentally, it can be 
argued [lit ] that the increasing particle energy entails a concomi- 
tant decrease in the ui-field energy, so total energy is conserved. 



A. Big Packet 

Consider a spread-out, real, positive, initial wave func- 
tion such as 

W^ = (2^5W e ^> (45) 
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where D >> a. We shall neglect the effect of the Hamil- 
tonian. We shall see that the wave function collapses 
fairly rapidly to an approximately spherical wave func- 
tion of size a, center location consistent with the Born 
Rule, and thereafter collapses more and more slowly to 
a smaller and smaller radius. 

First we calculate the ensemble average of the operator 
O = |X)(X|, so O is the ensemble probability density at 
X. For any density matrix, it follows from the collapse 
part of Eq.gOj), using Eq.(|5g|. that 



d|X)(X| 



dt 



-AiV 2 (X| 1 - e -l^[ Xi - x «] p(t)\X) = 0. 

(46) 

This, of course, doesn't say that collapse occurs but, if 
there is collapse, it says that the ensemble position prob- 
ability distribution does not change from the initial dis- 
tribution (Born Rule). 

In order to see that there is indeed collapse, consider 
the ensemble average of the modular momentum [l4l| op- 
erator O = cos P • nL, where n is a unit vector pointing in 
some direction. This is 1/2 the sum of two operators, one 
which translates the wave function by distance L in the 
n direction and the other in the — n direction. Thus, its 
expectation value gives the overlap of the wave function 
with itself (all the wave functions are real and positive) 
when translated. For any density matrix, it follows from 
the collapse part of Eq.(|40|) that 



cfcos P • nL 



dt 



= —\N' 



cos P • nL 



(47) 



Thus, for L >> a, the ensemble average of the overlap 
rate of the collapsing wave functions decreases as w AiV 2 
but then it slows, e.g., for L = a, the collapse rate is 
« .2XN 2 . 



B. 



Small Packet 



If the size of the wave function is less than a, one can 
utilize an approximate density matrix evolution equation 
obtained by expanding the exponential in Eq. (|40j) . retain- 
ing only the leading term: 



dp® 
dt 



2MN 



,/>(*) 



XN 2 
~4oJ 



Y^[Xi,[Xi,p(t)]]. (48) 



The state vector evolution which yields this density ma- 
trix evolution is 2 

|</,, i) = re^^'™^ 

. e -zx /o dt ' EtiK^O-VsAWa- 1 ^] 2 !^ o)(49) 

When the initial wave function is a gaussian, such as 
Eq. (|45|) with D < a, since the Schrodinger equation is 
quadratic in P and X, the solution is a gaussian. The 
exact solution to this problem can be found. Hi|[Ifi We 
shall arrive at it here using the formalism we have pre- 
sented. It suffices to solve the one-dimensional problem 
since, with initial wave function (|45l) , Eq. (14T)1) is the prod- 
uct of three terms, one for each dimension. 

We assume that the wave function at any time has the 
form 

^(X,t) = e- AWx2+SWX+C(t) , (50) 

and proceed to solve the Schrodinger equation which fol- 
lows from the time derivative of (|4"9l) : 



d_ 

di 



d 2 



2m dX 2 
X 



A 2 



4A 



w 2 (t) - T w(t)X + —X 2 ip(X, t), (51) 



A 



A 



where m = NM and A = XN/\/2a. Inserting (1501 into 
([5Tj) we obtain 



A = 
B = 



-2i 



-A 2 



-2i 



AB 



A 2 
A ' 
A 
" A 



w(t). 



(52a) 
(52b) 



Eq. (|52al) is a Ricatti equation, and can be solved by 
the ansatz A = (m/2i)F/F. It follows from (|52al) that 
F = F(2i\ 2 /m\). Thus, F = exp±a(l + i)t, where 
a = A/a/ mX and 

ma(l - i) e Q ( 1 +' i )* - K e - a( - 1+t ^ 

2 e *{i+i)t + Ke - a (i+i)t ■ ( - 53 - ) 

where K is a constant depending upon D. 

We see that A, which characterizes the squared stan- 
dard deviation of X, is the same for all w(t). Thus, the 
wave function approaches an equilibrium size, indepen- 
dently of its initial spread D. The equilibrium occurs 
because the Schrodinger evolution tends to spread the 



2 How can Eq. l|49j l. where w is just a function of t, arise from 
Eg. (I401 I. where w is a field, depending upon x as well as t? As 
far as I am aware, this has not been discussed before, so we 
treat it in Appendix [B] More generally, it involves changing the 
collapse-generating operators A a to a new, equivalent set, with 
concomitant change of white noise functions w a (t) to a new, 
equivalent set. 
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wave function while the collapse evolution tends to nar- 
row it. This takes place in characteristic time a -1 ~ 
5 x 10 4 /7Vsec. The equilibrium spread in X (its standard 
deviation) is l/y/2(A + A*) = l/V2ma « 4/A^ 1 / 2 cm. 

We shall henceforth assume either that the collapse 
process starts at negative times so that equilibrium is 
reached at time 0, or that A initially has its equilibrium 
value. Putting that value into (|52b[) gives 



B = -a(l + i)B + ^w(t), 

A 



(54) 



with solution 



B(t) = v / dt / w(t , )e-^ 1+i ^ t - t '\ (55) 
A Jo 

Knowing A and B, the expectation values of position and 
squared position can be found from (I5"U)) : 

^ (il>,t\X\1>,t) _ B + B* 



foM^M) 2(A + A*) 

2 / ^(tOe-^*-*') cos a(t-t'), (56a) 
A Jo 

1 



^ ' ~ (iMIM [ 1 + 2(A + A*) 



. (56b) 



To complete the solution, we need to find C(t) but, 
since it is used to find the probability density (^f),t\tp,t), 
it is best that we calculate that directly from the 
Schrodinger equation: 

-2^,t\X 2 \i>,t) 
1 [w{t)-2\{X)] 2 {^M,t) 



2A 



A 2 



1 



A A + A 



■(V,#,t). (57) 



Therefore, omitting the time-dependent factor arising 
from the last term of ff5T|) (which is absorbed in the nor- 
malization of the probability), and defining a new set of 
white noise functions 

v(t) = w(t) - 2X{X) 

= w(t)-2a f dt'w{t')e- a(t - t ' ) cos a (t -*')(£ 8) 
Jo 

the probability density is simply 

(fl>,t\i>,t) = e -^J'o** / * a (*'). (59) 



We note that Dw = Dv, since it follows from (|58l) that 
the Jacobian of the transformation from iu's to v's has 
l's on the diagonal and 0's above the diagonal. 



In order to use (|59j) . it is necessary to obtain the inverse 
of the transformation (|58"|) . This can be done by taking 
the second derivative of (l58l). with the result 



(60) 



<Pw(t) d 2 v(t) , 9 dv(t) 

+ 2a — - h 2av(t). 



dt 2 



dt 2 



(II 



Defining v(t)'s Brownian motion B(t) by v(t) = 
dB(t)/dt. It then follows from (|60| that w(t) can var- 
iously be written as 



w(t) = v(t)+2a / dtxv(tx) 
Jo 

dh 



+2a' 
ft 



n 



dt 2 v(t 2 ), (61a) 



o 



v(t) 



2a / dtiv{h)[l +a(t-ti)], (61b) 
Jo 

v(t) + 2aB{t) + 2a 2 [ dt'B(t'). (61c) 



It then follows from the first equation in ([SB"]) that (X) 
can be written as 



(X) 



1 



'mX 



B(t) 



dt'B(t') 



(62) 



One can then show, using ([50)1 . (|55|) . (|61aj) and (|52^) . that 



(P) = 2iA(X) - iB{t) = ^B(t). 

A 



(63) 



This problem is completely solved. We see from 
Eqs. (|6"2")) . (l()3")) that, after the equilibrium packet size is 
achieved, the momentum expectation value undergoes 
Brownian motion and the position expectation value un- 
dergoes a motion that can be described as Brownian-)-. 

Any expectation value can be calculated, and any en- 
semble average expectation value can be calculated. For 
example, although it can readily be found using the den- 
sity matrix, the ensemble average of the squared posi- 
tion expectation value can be found from Eq. (l56bp , us- 
ing (X) 



v(t)v(t') 



(l/2A)[u;(i) -«(<)] (Eq.flSB)), Eq.(j6Tb]), and 
X6(t - t') (which follows from (|59"]l): 



(X 2 ) 



2ma 


mX 


1 


^! 


2ma 


m J 


1 


m 


2ma 



dt'v(t')[l + a(t - t')] 

- [ dt'[l + a(t-t')] 2 
71 Jo 



,2+3 



f/3]. 



(64) 



X ~ t behavior occurs for classical Brownian motion, 
modeled as a particle undergoing Newtonian dynamics 
with a random force and a viscous damping force. In 
this case, the average Brownian 'step' size is constant in 
time. 
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X ~ t 3 behavior occurs for classical Brownian motion 
when the viscous damping is removed. This is essentially 
because the average Brownian 'step' size increases with 
time. 

So, we have the picture of the final result of collapse, 
a wave packet of equilibrium size which undergoes classi- 
cal random walk without viscous damping, with momen- 
tum generally increasing as it undergoes classical random 
walk. 



VII. COLLAPSE OF INTERFERING PACKETS 

It follows from the density matrix evolution Eq. (|4T))) 
that the interaction picture density matrix p(t) = 
W{t)p{t)U{t) (U(t) = exp(-iP 2 /2m)) satisfies 



dp(t) 



dt 

with solution 
p(t) =Te 



XN 2 U^(t) 



[Xi-X, 



U(t)p(t), (65) 



-\N 2 / Q 4 dt' U f (t') 



U(t') 



or, going back to the density matrix p(t), 



P(0) 
(66) 



p(t) = Te 



-\N 2 Jg dt'U(t-t') 



1— e 4a 



uHt-t') 



Po(t). 
(67) 

where po(t) = exp(— iHt)p(0) exp(iHt) is the density 
matrix without collapse. In the position representation, 
Eq.® is 



r z 2 («-«') 

-\N 2 f* dt' [l-e t^ 2 



<X|p(i)|X') = Te 



Z(t-t r ) = (X-^— ^V) 



fx' 



J (X|poW|X'), 
(68a) 

t —^V). (68b) 



We now note that, because [Xi — X[, Vj + Vl-] = 0, it 
follows that [Z(t - t') 7 Z(t - if')] = 0, and so the time- 



ordering operation T may be removed from Eq. (|68a| ), 
Also because this commutator vanishes, any product of 
powers of Z's can be written in 'normal-ordered form,' 
by which we mean that the X's are to the left of the 
Vs. Denoting the normal ordered form by : :, Eq. (|68a[) 
becomes 



Z 2 (t-t') 

l-e 4^ 



(X|po(t)|X') 
(69) 

We shall apply Eq. dMl) to the case where the uncol- 
lapsed density matrix po(t) is constructed from a number 
of wave packets, 

(X\p (t)\X') =J2cnCct>n(X,t)cj>* n ,(X',t). (70) 



The wave packets cp n (K,t) are to have well-defined mo- 
menta k„(X) at (almost) each point of the wave packet, 



which itself has dimensions large compared to the wave- 
length. Thus, a wave packet could be a laboratory 'plane 
wave,' a good approximation to an eigenstate of momen- 
tum k. It could be a cylindrical wave packet or a spheri- 
cal wave packet of momentum magnitude k such as might 
be obtained by putting the 'plane' wave packet through 
a slit or a circular hole. 

An important feature of such a packet <j>j (X, t) is that 



(X|P|0„,t) = -V0 n (X,i) » k„(X)^(X,t) (71) 

is a very good approximation. Another important feature 
of such a packet is that (almost) each point in each wave 
packet can be considered as moving on a straight-line 
trajectory with constant velocity k„(X)/m. 

Putting together Eqs.([6l ^ .(|6 % ([70 ]) . we obtain for the 
ensemble's probability density at X: 

(X\p(t)\X) = £ c „<£,0 n (X,t)&,(X,t) 



So dt' [ 



\N 2 J*dt' l-e 4a 



X n (t-t')-X z(t-t') 



where 



x n (t-0^x-^ffi(i-t'). 



(72) 



(73) 



That is, consider a point on the nth packet which is lo- 
cated at X at time t. Then, X n [t — if) is the location 
that point had on the nth packet at the earlier time if . 

To summarize, we have seen in Eq. pp|) or (|44p that, 
when a clump is put into a superposition of two places 
with constant separation D, the two states play the gam- 
bler's ruin game, so that the off-diagonal elements of the 
density matrix decay at the rate A-/V 2 [l — cxp — (D 2 /4a 2 )]. 
Eq. (l75)l says that, for a superposition of packets, the 
points on the packets, which end up at the same place 
X at time t, may be thought of as playing the gambler's 
ruin game with each other on the way to X, with the 
above-mentioned distance-determining rate now varying 
with time, governing the collapse all along the way. 

Although it is not of concern here, we mention that, of 
course, the spatially separated points of a single packet, 
or of multiple packets, likewise mutually play the gam- 
bler's ruin game, and that description is obtained by con- 
sidering the off-diagonal elements of the density matrix. 



A. Mach-Zender Interference 

As is well known, the Mach-Zender interferometer has 
a rectangular shape, say, with half-silvered beam-splitters 
at the lower left and upper right corners, and fully- 
silvered mirrors at the other two corners. An incoming 
wave packet splits into two equal packets at the first beam 
splitter. The packet going — >, j" has its sign reversed 
when it reflects at 90° from the front-surfaced mirror. 
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The packet going f, — h t gets no ne t sig n change: one 
sign change at the first, front-surfaced, beam splitter, one 
at the front-surfaced mirror, none at the second, back- 
surfaced, beam splitter. Thus, without collapse, there is 
no output in the j" direction. 

Although there are certainly velocity changes of the 
packets, they take place over a relatively brief time in- 
terval, so Eq.(|72|) may be applied seriatim. Let t be the 
time interval separating emergence from the two beam 
splitters. At time 0, the two packets start off with 
ci = c 2 = 1/V2. Thereafter, \X 1 {t-t'yX 2 (t-t')\ » a. 
Moreover, the collapse rate is unaffected if a packet 
changes sign. Finally, at time t, the second beam split- 
ter has just made the amplitudes Ci = c-i = 1/2 and, if 
packet 1 is 0i(X, t), then packet 2 is — </>i(X, t). 

Therefore, in Eq.([?2j), since exp -[X„ - X„] 2 /4a 2 = 1, 
cxp — [Xi — X2] 2 /4a 2 « 0, for a point X in the superposed 
wave packets, 



(X|p(t)|X) = ^ 1 (X,t)| 2 2[l-e 



(74) 



Eq.([72D yields 



(X\p(t)\X) = ^-[2 + (e 2Me 



-\N 2 [t~ f* dt' 



-2ikb6\ 



A 



1 + cos(2kbd)e 



-MN 2 [l-f f ° dye-" ] 



2A 2 cos 2 (kb6)e 



XtN 2 [l-^ X {b/a)] 

-\tN 2 [l-^ X (b/a)] 



-A 



1 



(76) 



where %(6/a) = erf(6/a). 

Thus, we see that the two-slit two-packet interference 
pattern decays while the single packet non-interference 
pattern builds up as time increases. For b >> a, the 
packet separation is >> a for almost all the time and 
the collapse rate is AiV 2 , as in the previous section. For 
b « a, the collapse rate is XN 2 b 2 /3a 2 . 

This concludes our discussion of free particle collapse 
dynamics. 



and the probability that the clump emerges in the upward 
direction is 



/> = / rfX<X|p(t)|X) = i[l-e-^ 2 *]. (75) 



We see that, as time spent in the interferometer increases, 
P-f — s- 1/2 since, asymptotically, only one packet survives 
to hit the second beam splitter, and that packet has equal 
likelihood of going | or — >. 



B. Two Slit Interference 

It should be clear that Eq. ([72l can be applied to any 
interference or diffraction problem. Here we shall just 
consider the effect of collapse on the two-slit Fraunhofcr 
interference pattern, with neglect of single-slit diffrac- 
tion. 

Take the two slits to be located at x = ±6 and parallel 
to the z-axis. The pattern is observed at a point X on a 
screen located at y = L >> b. The vectors from the slits 
to the point on the screen are 1*12 = X ± ib. Defining 8 
as the angle X makes with the y-axis, and working only 
to first order in 0, then X « and n,2 ~ |X|±60. 

The two cylindrical packets </>i.2 have wave number k, 
and are of equal amplitude ci ; 2 = l/v2 when they em- 
anate from the two slits. They reach the screen at point 
X at time t, where their amplitudes are Aexp±ikb9. 
We are to consider the parts of the two packets trav- 
eling with speed k/m on a straight line, each from its 
slit to the screen. Their separation at any time t' is 
\Xi(t-t')-X. 2 (t-t')\ =2b[l-(t'/t)\. Putting this into 



Appendix A: Proof That R and S must be diagonal 

We prove here that the real symmetric operators R 
and S in the Stratonovich Schrodinger equation for the 
un-normalized state vector, 



,t) = [MB' + Sdt]\4>,t) 



(Al) 



must be diagonal in the \a n ) basis. This is in order 
that Eq. (|Al[) give rise to the ltd gambler's ruin condi- 
tion Eq.©, 

dx n (t) = b n (x)dB(t). (A2) 

After putting Eq.flU), dB' = dB + fdt, into Eq. (|ATl) . 
we convert that Stratonovich equation to an Ito equation, 
with the result 

d\((>, t) = [RdB + Vdt]\<f>, t) where V = S + Rf + ^R 2 . 

(A3) 

We note that V is also a real symmetric operator and, if 
we show R and V must be diagonal, then S must also be 
diagonal. 

Using the rules for manipulating Ito equations, it is 
straightforward to find 

d|<M)(<M = {[RdB + Vdt],\4>,t){<f>,t\} 

+\dtR\<f>,t){<f>,t\R, (A4a) 

d{<f>, t\<p, t) = 2[{<(>, t\R\4>, t)dB + {(j>, t\V\<j>, t)dt] 

+\dt(<j),t\R 2 \(l) 1 t)((l) 1 t\R, (A4b) 

where {M, N} = MN + NA1. Defining the density ma- 
trix p(t) = \<j),t){<t>,t\/{(t>,t\<t>,t) and M = TraceMp, we 
obtain from Eqs.© and the Ito rules: 

dp = [{R, p} - 2pR]dB + dt[{V, p} - 2pV] + 

Xdt \[RpR - pW] - 2R[{R, p} ~ 2pR}} . (A5) 
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Now, x n {t) = (a n \p(t)\a n ) . Thus, in order that the 
diagonal elements of Eq. (IA5|) agree with Eq. (|A2[) , we see 
that the diagonal elements of Eq. (|A5j) which do not mul- 
tiply dB must vanish for arbitrary p: 

= [{V, P } nn - 2 Pnn V] + X [[(RpR) nn - Pnn W] 

-2R[{R, p} nn - 2p„„i?]] . (A6) 

where M nrn = (a n \M\a m ) 

First, suppose that p mm — 1, where m ^ n, and 
all other matrix elements of p vanish. It follows from 
Eq. ([A"6f that 



= (RpR) r 



(R„ 



(A7) 



That is, all the off-diagonal elements of R vanish, so R is 
diagonal. 

Second, choose a density matrix for which p nn , p mm = 
1 — Pnn, Pnm do not vanish, but all other matrix elements 
of p do vanish. Then, using the diagonal nature of R, 
Eq. (|A6p may be written as 

2V nm p nrn \\ 2p n7l \ 

~t~Pnn[l Pnn] 2iV nn V mrn ) 
-|-A(i? rm Rmm)[(Rnn ~l~ Rmm) 

^{RnnPnn + Rmrn (1 — pnn)) ■ (A8) 

For fixed p nn , a viable density matrix (non- negative 
eigenvalues which add up to 1) exists for \p nm \ < 
\J Pnn{\ — Pnn)- But, as p nm is varied, the first term in 
Eq. (IA8l) varies while the rest of the terms remain fixed. 
Thus, the first term must vanish, and this means that 
V nm = for n 7^ m, i.e., V is diagonal as well as R. 



Appendix B: Transformation of Operators and 
White Noise 

Consider the general CSL form for the evolution of the 
state vector, Eq. (l33|) 

\^ t )=Te-fo dt 'i iH ( t ')+ix^J^')-^^}^ (B1) 

We introduce a real orthonormal set of vectors u'g, i.e., 

Sq u u /3' = <W'> Y,i3 u< fi u i = <W- Defining white 
noise functions v^(t) and complete commuting set of 
operators Z 13 by w a {t) = Y^^ u ^{t) and A a (t) = 
'Yl,pU c pZP(t), one readily sees that, in the exponent of 
Eq.pqi. 

^K(t) - 2XA a ] 2 = 5^[i^(t) - 2XZ f) } 2 . (B2) 



The Jacobian of the transformation from id's to v's is 1 
so, in using the Probability Rule (|32b|) . Dw = Dv. 



We wish to apply such a transformation to Eqs. (l4TT) . 
()42|) which, for simplicity, we limit to one-dimensional 
space: 



\<j>,t) =Te- l K dt '^w p2 

. e ~ 4T /o dt' fdx'[w(x',t')~2XNA(x')} 2 



A(x') 



1 



(tto 2 ) 1 / 



-e 2 



,0),(B3a) 
(B3b) 



We shall use as orthonormal functions the harmonic os- 
cillator wave functions 

1 2 1 

u n {x) = C n H n ix/a)e ^ x where C n = =. 

Vir 1 / 2 2 n n\a 
(B4) 

With the definitions u n (i) = J dxw(x,t)u n ix) and Z n 
| dxAix)u n ix), the exponent in Eq. (|B3a[) may be written 
as 



1 

4A 



dt' J dx'[w{x',if) - 2XNAix')] 2 = 
-ttV / dt'[v n it')~2XNZ n } 2 . (B5) 



Thus, we see how a white-noise field gets converted to an 
equivalent sum of white noise functions. 

Using the identity exp(— t 2 + 2tz) = X)^Lo t n H n iz)/n\, 
with t = X/2a, z = x' /a, we find 



dx 



(Tra 2 ) 1 /* 



e 2 



C n (7ra 2 )V4 



m=0 



U„(l) 

(A > /2a) r 



e 



?2 (A > / V / 2a) r 



dxu m (a;)u„(a;) 

(B6) 



This leads to the density matrix evolution equation 



dpjt) 
dt 



p2 



-.Pit) 



XN 2 



2MN 
_x* iX/V2a 



71=0 



(B7) 



If we expand exp— X 2 /4a 2 , we see that the n = term 
goes as iX/a) 4 and the rest of the terms go as (X/a)™ 
to lowest order. Therefore, the lowest order term comes 
from n — 1. Upon neglect of the higher order terms, this 
gives the density matrix evolution equation 



dpjt) 
dt 



P 2 



2MN 



XN 2 



[X,[X,p(t)]]. (B8) 



which is identical to the one-dimensional version of 
Eqj 
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